Entanglement dynamics of three-qubit states in many-sided noisy channels 
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We study entanglement dynamics of pure three-qubit Greenberger-Horne-Zeilinger-type (GHZ- 
type) entangled states when one, two or three qubits being subjected to general local noise. Em- 
ploying a lower bound for three-qubit concurrence as an entanglement measure, we show that for 
some many-sided noisy channels the entanglement dynamics can be completely described by the 
evolution of the entangled states in single-sided channels. 
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I. INTRODUCTION 

It is widely accepted nowadays that entangled states 
of multiparticle systems are the most promising resource 
for quantum information processing. At the same time, 
entanglement of complex systems is known to be very 
fragile with regard to decoherence, which may appear, 
for instance, due to a transmission of the whole quantum 
system or some its subsystems through communication 
channels. In practice, moreover, it is often required to 
distribute parts of an entangled multipartite system be- 
tween several remote recipients. In this case, each sub- 
system is coupled locally with its environment. Such a 
coupling of a quantum subsystem to some environmen- 
tal channel leads to decoherence of the multipartite sys- 
tem and usually to some loss of entanglement. For suc- 
cessful practical utilization of entangled states, it is of 
great importance to quantify entanglement of the com- 
plex quantum systems and describe quantitatively their 
entanglement dynamics under the action of local (noisy) 
channels. 

There are two main approaches in the literature to in- 
vestigate entanglement dynamics of multiparticle entan- 
gled states. First one is to study state evolution of par- 
ticular (usually maximally entangled) states under the 
action of some chosen noisy channels and deduce entan- 
glement dynamics from the state evolution [H-Q. This 
approach, however, is restricted by our choice of the en- 
tangled states and the noise models. Recently, a com- 
pletely different approach for the description of entan- 
glement dynamics, which is based on evolution equation 
for entanglement, has been developed This latter 

approach allows to obtain a direct relationship between 
the initial and the final entanglement of the system in 
which just one of its subsystems is subjected to an arbi- 
trary noise. Unfortunately, the suggested evolution equa- 
tions can not be straightforwardly generalized to the case 
when more than one subsystem undergo the action of lo- 
cal noisy channels. 

In this work, we investigate entanglement dynamics of 



initially pure three-qubit states when one, two or three 
qubits are affected by general local singlc-qubit noisy 
channels. At first, we shall consider the entanglement dy- 
namics of the maximally entangled Grccnbcrgcr-Hornc- 
Zcilingcr (GHZ) state, which can be written, in the com- 
putational basis, as 



|GHZ> = i= (|000) + |111» 



(1) 



Later, we shall generalize our discussion to all (GHZ- 
type) pure three-qubit states which can be obtained from 
the state ([T]) by local unitary transformations. To de- 
scribe the influence of the local noisy channels on the en- 
tangled states we shall use quantum operation formalism 
Since there is no an analytically computable mea- 
sure of entanglement for multiqubit states [lfj, a lower 
bound for multiqubit concurrence [ll[ will be utilized to 
access the entanglement dynamics. We shall show that 
for some noisy channels the complex two- and three-sided 
entanglement dynamics, i.e. when two or three qubits 
are affected by local noise, can be completely described 
by the evolution of the entangled system in single-sided 
channels (when just one qubit is subjected to local noise). 

This work is organized as follows. In the next section, 
we present the entanglement measure of use, the lower 
bound for multiqubit concurrence, and recall its proper- 
ties. In Sec. IIII1 we step by step analyze cases when 
one, two or three qubits of the three-qubit system are 
subjected to local noisy channels. A summary is drawn 
in Sec. HV1 



II. 



THE ENTANGLEMENT MEASURE 
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Concurrence, originally suggested by Wootters [12| to 
describe entanglement of an arbitrary state of two qubits, 
has been recognized as a very powerful measure of entan- 
glement. Although various extensions of the concurrence 
to the case of bipartite states, if the dimensions of the 
associated Hilbert (sub-) spaces are larger than two, has 
been suggested [L3l - [l5l |. a full generalization of the con- 
currence towards multiparticle states still remains chal- 
lenging [Io| . However, a formal extension of the concur- 
rence to multipartite case can be successfully approxi- 
mated by an analytically computable function, so-called 
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lower bound, which never exceeds such multipartite con- 
currence, but nonetheless is close enough to its values. 
To date, several lower bounds for the multipartite con- 
currence have been proposed [l|, [ltj [H| • Let us recall and 
exploit the lower bound for the multiqubit concurrence 
as suggested by Li et al. . 

The multiqubit concurrence for a pure three-qubit 
state \ij)) is given by 



c 3 (m = 



\ 



(2) 



where the pt = Tr \tp) denote the reduced density ma- 
trix of the i-th qubit which is obtained by tracing out 
the remaining two qubits. The concurrence for an arbi- 
trary mixed thrcc-qubit state can be formally defined by 
means of the so-called convex roof 



C 3 (p)= min Jjp, CM)) 



(3) 



relying on the fact that any mixed state can be ex- 
pressed as a convex sum of some pure states 
p = J^i Pi IV'i) (V'il- However, the minimum in this ex- 
pression is to be found among all possible decompositions 
of p into pure states \ipi). No solution has been found so 
far to optimize the concurrence © analytically [Io| . 

A simple analytically computable lower bound r^{p) 
for three-qubit concurrence can be given in terms of the 
three bipartite concurrences C ab l c (a,b,c = 1,2,3 and 
a^6^c^a) as 



\ 



E 

k=l 



(c ,l 2 |3 )2 + (c ,13|2 )2 + (c ,23|l )2 _ (4) 



Here, each bipartite concurrence C ab ' c corresponds to a 
possible (bipartite) cut of the three-qubit system in which 
just one of the qubits is discriminated from the other two 
qubits. For a separation ab\c, the bipartite concurrence 
(jab\c j g gj ven by a sum of six terms Ck which are ex- 
pressed as 

Cf |c = max(0,A^-A2-Af.-A 4 fe ) , (5) 



and where the A 



k , m, — 1..4 are the square roots of 

~ able 



the four nonvanishing eigenvalues of the matrix pp k , 

if taken in decreasing order. These matrices p p^ 
are formed by means of the density matrix p and its 
complex conjugate p* , and are further transformed by 



the operators {£ 



ab\c 
k 



- ab\c 



Lq, k = 1, 6} as: 



Pk = 5fc h ' C p*<S'fc h ' C - In this notation, moreover, Lq is the 
(single) generator of the group SO(2) which is given by 
the second Pauli matrix a y = —i (|0) (1| + |1) (0|); while 

the {i^ b ' c } are the six generators of the group SO(4) 
that can be expressed explicitly by means of the totally 
antisymmetric Lcvi-Ccvita symbol in four dimensions as 
{Lkl)mn = -iSkimn, k,l,m,n = 1..4 [HI. 



Since the lower bound ([4]) is just an approximation for 
the convex roof ([3]), it is of great importance to know 
how accurate this bound is with respect to the convex 
roof for multiqubit concurrence. It has been checked [TtJ 
by sampling 100 randomly generated density matrices p 
that the lower bound T3 (p) coincides with the numerically 
simulated convex roof C3 (p) for all density matrices with 
rank r < 4. This result has found its explanation in 
the fact that the lower bound r 3 (p), by its construction 
[ll|, relies only on four nonvanishing eigenvalues of the 

matrix pp^ h ' c , while this matrix may have at most eight 
eigenvalues. 

Taking into account the mentioned property of the 
lower bound T3(p), it is convenient to define a function 
of at most four variables 

f(w,x,y,z) = 

min (0, 2 max (w, x,y,z) — w — x — y~z) . (6) 



It is easy to see that C 



ab\c 



m,Xl,4,4), where 



A™, m = 1..4 are the square roots of the eigenvalues of 
the matrix pp^ c . Whenever the matrix pp^ fc ' c has less 
than four, e.g. two, eigenvalues, we shall use the function 
of just two inputs f(w,x) = f(w, x, 0, 0). 



III. ENTANGLEMENT DYNAMICS IN LOCAL 
NOISY CHANNELS 

Quantum operation formalism is a very general and 
prominent tool to describe how a quantum system has 
been influenced by its environment. According to this 
formalism the final state of the quantum system, that is 
coupled to some environmental channel, can be obtained 
from its initial state with the help of (Kraus) operators 



Pfin 



E 



Ki p ini K. 



(7) 



and the condition K\ Ki = I is fulfilled. Note that 
we consider only such system-environment interactions 
that can be associated with completely positive trace- 
preserving maps 

If the quantum system of interest consists of just single 
qubit which is subjected to some environmental channel 
A, then an arbitrary quantum operation can be expressed 
with the help of at most four operators @ . Let us define 
the four operators through the Pauli matrices as 



ai | 




! 


V2 \ 










—i 


V2\ 








V2\ 



where are real parameters and the condition 
Eli o? = 1 holds. 



3 



For threc-qubit system, local interactions of the qubits 
with channels A, B and C can be described by opera- 
tors, which arc constructed as tensor products of the 
single-qubit operators Ki(ai), Ki(bi) and Ki(ci). There- 
fore, if the three qubits are affected by local noise simul- 
taneously, the final state of the system ([7]) can be ob- 
tained from its initial state with the help of 64 operators 
Ki(a,i)<8>Kj(bj)<3)Ki(ci) = 1...4). If just two qubits 

undergo the action of local channels, the total number of 
Kraus operators is 16. 

A. Single-sided channels 

Suppose, just one qubit of the three qubit system, 
which is initially prepared in pure GHZ state ([T]), is sub- 
jected to a noisy channel A. The final state of the system 
is in general mixed and can be expressed by a rank-4 den- 
sity matrix p. This matrix is obtained from Eq. ([7} using 
the four operators 1® l®Ki{ai). Here we assumed, with- 
out loss of generality, that the third qubit undergoes the 
action of the channel. From the matrix p we can com- 
pute the three bipartite concurrences C 12 ' 3 , C 13 ' 2 , C 23 ' 1 
and the lower bound Q. The concurrences are given by 

C 12 ' 3 = /(a 2 , a 2 , a 2 , a 2 ), (9) 

c 13 ' 2 = y7 2 (a 2 ,« 2 ) + / 2 (4«I.), (io) 

C 23 ' 1 = ^p(alal) + p(al4), (11) 

where we dropped the arguments of the concurrences, 
e.g. C 12 ' 3 = C 12 l 3 ([l <g> 1 <g> A] \GHZ) (GHZ\) 

Interestingly, if the channel A is the bit-flip (ai ^ 
0,a2 7^ 0,03 = a4 = 0) or the bit-phase-flip (a\ ^ 
0,a 3 ^ 0,a 2 = a 4 = 0) channel [|, the entanglement of 
the three-qubit system, described by means of the lower 
bound T3(/o), never vanish. 



The description of the entanglement dynamics of the 
GHZ state (JT}) in single-sided channels can be general- 
ized to the case of an arbitrary pure threc-qubit state \ip) 
which can be obtained from the GHZ state by local uni- 
tary transformations. Such an extension is possible due 
to recently suggested evolution equation for the lower 
bound T3 [8J, which manifests that the entanglement dy- 
namics of an arbitrary pure state \ip) of a three-qubit 
system, when one of its qubits undergoes the action of 
an arbitrary noisy channel A, is subordinated to the dy- 
namics of the maximally entangled state, i.e. 

r 3 [(l®l® A)\il>)(il>\] = 

r 3 [(l <8> 1 ® A) \GHZ) (GHZ\] r 3 [|V>] ■ (12) 

Conclusively, the entanglement dynamics of an arbitrary 
pure three-qubit state, which is locally equivalent to the 
GHZ state and is subjected to an arbitrary single-sided 
channel A, is given by the lower bound Q, which is 
defined through Eqs. (|9"))- ([Tl"j) . and the evolution equation 
(I18|) . Moreover, based on the numerical simulations in 
[17j , we argue that the description of this particular case 
of the entanglement dynamics with the lower bound T3 
is as general as with the convex roof ([3]) for multiqubit 
concurrence. 



B. Two-sided channels 

If the three-qubit system is prepared into the GHZ 
state ((T|) and just two (let say, the second and the third) 
qubits are affected by local channels A and B simultane- 
ously, the final state of the system is given by a rank-8 
density matrix p. In this case, the three bipartite con- 
currences can be expressed as 



C 12 ' 3 = ^P(a\b\ + ajbl a\b\ + a§6§, a\b\ + a|fe|, a\b\ + a|&|) + p{a\b\ + a?6§, a\b\ + a\b\, a\b\ + a\b\, a\b\ + a\b\) , (13) 
C 13 ' 2 = ^P{a\b\ + afbl albl + a\b\, a\b\ + a§6§, a\b\ + a\b\) + f\a\b\ + a\b\, a\b\ + a\b\, a\b\ + a\b\, a\b\ + a\b\) , (14) 
C* 2311 = ^P{a\b\ + a\b\, a\b\ + a\b\, a\b\ + a\b\, a\b\ + a\bl) + f 2 (4b 2 2 + a\b\, a\b\ + a§b§, a\b\ + a\b% a\b\ + a\b^) . (15) 



Although no general conclusion about the entanglement 
dynamics of the three-qubit system can be made from 
the structure of these bipartite concurrences, there are 
some interesting particular cases. If each of the channels 
A and B is just the bit-flip or the bit-phase-flip channel, 



the squared lower bound r| is given by 
i[/ 2 ( a 2 ,a 2 ) + / 2 (6 2 ,& 2 ) + / 2 (a 2 ,a 2 )/ 2 (& 2 ,6 2 )] , (16) 

where i,j = 2,3. Symbolically Eq. (|16j) can be written 

as 

tI{[1® A®B]p) = \ (A 2 + B 2 + A 2 B 2 ) , (17) 
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where A = C 12 l 3 ([l ® 1 <g> A] p), B = C 12 l 3 ([l ® 1 ® 5] 
and /o = \GHZ) (GHZ\. Other words, entanglement dy- 
namics of the GHZ state in two-sided noisy channels, 
when the channels A and B are the bit-flip or the bit- 
phase- flip channels, is completely defined by the dynam- 
ics of the state in the single-sided channels. Although this 
result holds only for mentioned channels, it can be ex- 
tended to the case of an arbitrary pure three-qubit state 
p = \ip) (-01 which is locally equivalent to the GHZ state. 
This extension is based on the evolution equation for bi- 
partite concurrence @ which has very similar sense and 
structure to the evolution equation (TT8"]) for the lower 
bound. For the bipartite concurrence C 12 ' 3 , for example, 
the evolution equation is given by 

C 12|3 [(1®1® A)\^){^\] = 

C 12|3 [(l g> 1 g> A) \GHZ) (GHZ\] C 12|3 [|V>] . (18) 

It is important to note that Eqs. (|13 p -(|15 |l are com- 
puted for a rank-8 density matrix and, therefore, the 
lower bound T3 obtained with the help of these equations 
may differ from the convex roof ([3]) for the concurrence 
significantly. However, due to the assumption that the 
channels A and B arc the bit-flip or the bit-phase-flip 
channels, the state [1 (g> A ® B] \GHZ) (GHZ\ has rank 
four. Therefore, it can be argued that Eq. (|T7|) . formu- 
lated for the lower bound T3, remains valid if the convex 
roof ([3]) is substituted in its left hand side (lhs). 

C. Three-sided channels 

In the previous subsection we have analyzed the case 
when just two qubits of the three-qubit system are sub- 
jected to local noisy channels. Similar analysis can be 
made when all three qubits are affected by local chan- 
nels A, B and C simultaneously. If the system is initially 
prepared in the GHZ state (J]), the final state density 
matrix, derived from Eq. (J7J), has rank eight. The ana- 
lytically computed concurrences C 12 ' 3 , C 13 ' 2 and C 23 ' 1 in 
this case have very complex structure and, therefore, we 
do not display them here. There are significant simplifi- 
cations in the structure of the concurrences and the lower 
bound if all A, B and C are bit-flip channels or if two of 
these channels arc bit-phasc-flip and the remaining one 
is a bit-flip. In these two cases, the squared lower bound 
r| can be written as 

tI{[A®B®C\ P ) = \ (A 2 B 2 +B 2 C 2 +A 2 C 2 ) , (19) 



where A = C 12 l 3 ([l ® 1 ® A] p), B = 
C* 12 l 3 ([l® l®B]p), C = C 12 l 3 ([l <g) 1® C]p) and 
p = \GHZ) (GHZ\. The entanglement dynamics of the 
GHZ state in the three-sided noisy channels is given just 
by the dynamics of the state in the single-sided channels. 
As in the previous section, this result can be generalized 
to the case of an arbitrary three-qubit state (which is 
locally equivalent to the GHZ state) due to the evolution 
equation for bipartite concurrence @. Moreover, the 
fact that the state [A®B®C] \GHZ) (GHZ\ has rank 
four, for mentioned channels, suggests to substitute the 
convex roof © in the lhs of Eq. (|19p without loss of any 
information about the entanglement dynamics. 



IV. SUMMARY 



Our analysis in the previous section suggests that 
the entanglement dynamics of a pure three-qubit GHZ- 
type entangled state \ip) in some local many-sided noisy 
channels can be completely described by evolution of 
the entangled state in single-sided channels. More 
specifically, if just two qubits are subjected to local 
channels A and B, so that the channel [1 <g> A <g) B] 
is a combination of just bit-flip and/or bit-phasc- 
flip channels, the lower bound for three-qubit concur- 
rence r 3 ([l <g> A (g) B] \ip) (ip\) can be expressed without 
loss of generality through the bipartite concurrences 
C 12 l 3 ([l ® 1 ® A] \ip) (tp\) and C 12 l 3 ([l® l®B]\if>){il>\) 
by Eq. (p~T|) . If all three qubits undergo the action of lo- 
cal noisy channels A, B and C, which are neither bit-flip 
channels or two of these channels are the bit-phase-flip 
and the remaining one is the bit-flip, the lower bound 
Tz([A <g> B CS> C] \ip) (ip\) is again given just by the bipar- 
tite concurrences as it is displayed in Eq. (|19j) . 

It is important to note that there are only two nonva- 
nishing eigenvalues in the expression (j9]) for the bipartite 
concurrence C 12 ' 3 ([l ® 1 <8 A] \tp) (ip\), if the channel A is 
the bit-flip or the bit-phase-flip channel. Such a bipartite 
concurrence can be directly measured (3- Eqs. (TT7)) and 
(fT9"| allow one to access the complex many-sided entan- 
glement dynamics just by measuring the bipartite con- 
currence in the simplest case, when just one qubit is af- 
fected by a local noisy channel. 
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